Linear optical absorption spectra of mesoscopic structures in intense THz fields: 

free particle properties 
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We theoretically study the effect of THz radiation on the linear optical absorption spectra of semi- 
conductor structures. A general theoretical framework, based on non-equilibrium Green functions, 
is formulated, and applied to the calculation of linear optical absorption spectrum for several non- 
equilibrium mesoscopic structures. We show that a blue-shift occurs and sidebands appear in 
bulk-like structures, i.e., the dynamical Franz-Keldysh effect [A. -P. Jauho and K. Johnsen, Phys. 
Rev. Lett. 76, 4576 (1996)]. An analytic calculation leads to the prediction that in the case of 
superlattices distinct stable steps appear in the absorption spectrum when conditions for dynamical 
localization are met. 
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I. INTRODUCTION 

Light absorption can be described in terms of a process 
where polarization is induced in the medium. To linear 
order in the electric field component of the traversing 
light, E, the induced polarizability, V, can be expressed 
in terms of the dielectric susceptibility x ^ 



r(t) 



dt' x{t,t')£{t') 



If the absorbing medium is in a stationary state, the sus- 
ceptibility depends only on the difference of its time ar- 
guments, i.e., i') = x(t — t'). Under these condi- 
tions Maxwell's equation for £ is an algebraic equation 
in frequency space and one finds that the absorption is 
proportional to the imaginary part of x{^)- However, 
under non-equilibrium conditions, which are the topic of 
the present study, the susceptibility is a two-time func- 
tion, and Maxwell's equation remains an integral equa- 
tion even in the frequency domain. Further progress 
hinges upon two steps: first, one has to develop meth- 
ods to calculate the non-equilibrium susceptibility func- 
tion and second, one has to specify what sort of light 
wave or pulse is used in the absorption experiment. The 
present work addresses both of these problems. As far 
as the time-dependence of the probe pulse is concerned, 
two specific situations are examined. First, consider an 
undoped semiconductor placed in an intense THz field; 
we assume that the THz field is not able to induce po- 
larization, i.e., no carriers are excited in the conduction 
band. Such a system is in a non-equilibrium state, i.e., 
the susceptibility is a two time function. Properties of 
such systems have been investigated experimentally us- 
ing the Free Electron Laser (FEL) as a source for intense 
THz fields | ; many interesting properties have been dis- 
covered, and others predicted, such as photon-assisted 



tunneling ||^, dynamical localization and absolute nega- 
tive conductivity , ac Stark effect , dynamical Franz- 
Keldysh effect and formation of sidebands A sec- 
ond example consists of ultra-fast transients. Consider an 
undoped semiconductor structure subject to an external 
intense static field. At some time instant a population 
of carriers is pumped into the conduction band. These 
mobile charges will rearrange themselves so as to screen 
the external field. While the screening is building up 
the susceptibility of the system is a two time function. 
Using femtosecond laser techniques, which are able to 
probe the time scales in which screening is building up, 
experiments investigating the non-equilibrium properties 
of such systems have been performed 10 1. 

Bandgap engineering techniques of semiconductor 
compounds, such as Molecular Beam Epitaxy (MBE), 
allow spatial modulation of the bandgap down to atomic 
resolution. It is possible to break the translational sym- 
metries of bulk crystals, induce new ones and reduce the 
degrees of freedom with these techniques. Construction 
of systems which are effectively two dimensional (2D) and 
even one dimensional (ID) with regard to electron mobil- 
ity and optical properties is today a standard procedure. 
Another example of such man-made structures are super- 
lattices (SL), i.e., an engineered periodic potential in the 
growth direction of the sample. The interplay between 
the mesoscopic properties and the dynamical properties 
can lead to many interesting phenomena such as absolute 
negative resistance for the transport properties [ pd| and 
the rich features in the optical properties which are the 
subject of the present work. 

The purpose of this paper is to present a theoretical 
study of light absorption in mesoscopic systems subject 
to intense THz [far infrared (FIR)] fields. We consider 
undoped systems, which implies that there are no car- 
riers in the conduction band. Thus, the near infrared 
(NIR) inter-band absorption is the dominant absorption 
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process. The non-equilibrium nature of the system ne- 
cessitates the use of special theoretical tools; we have 
chosen to apply the non-equilibrium Green function tech- 
niques [ |l2yi3|| . In particular, this method allows us to 
treat the intense FIR field non-perturbatively, and de- 
fines a framework in which screening can be treated sys- 
tematically. Our analysis consists of the following steps. 
Starting from a two-band Hamiltonian we derive a for- 
mal expression for the inter-band susceptibility in terms 
of non-equilibrium Green functions. Next, we use the 
general expression to derive the NIR absorption spec- 
trum for non-interacting particles (Coulomb interactions 
will be discussed in a subsequent paper, see also below). 
Finally, we give explicit results for a number of special 
cases (3-, 2-, 1-d systems and super lattices), and discuss 
the physical implications. 

This paper is organized as follows. In section |l| we de- 
rive a general expression for the two-time dielectric inter- 
band susceptibility. Section III relates the susceptibility 
to the measured absorption, considering both continuous 
wave measurements as well as short white light pulses. 
The single-particle Green functions and the correspond- 
ing spectral functions, which determine the susceptibility, 
are defined in section IV and related to the generalized 
density of states 1^ (GDOS), which, in turnj^ is shown 
to determine the optical absorption. Section ^ considers 
bulk-like systems, and we obtain analytic results for the 
GDOS, which is analyzed in some detail in terms of the 
side-band picture. The properties of light absorption in 
superlattices are treated in section VI. Specific atten- 
tion is paid to conditions where dynamical localization 
occurs, and we show how it a ffects the absorption spec- 
trum. Finally, in section VIl we have some concluding 
remarks. 



II. THE DIELECTRIC INTER-BAND 
SUSCEPTIBILITY 

We shall now derive an expression for the dielectric 
inter-band susceptibility using non-equilibrium Green 
functions. The microscopic operator describing inter- 
band polarization is 



F(t) = ^4 [aiit)bk{t) 



bl{t)ak{t) 



(1) 



Here dk is the dipole matrix-element, a|,(t) [afc(i)] are 
the conduction band electron creation [annihilation] op- 
erators and bl,{t) [bk{t)] are the valence band creation 
[annihilation] operators. The linearized Hamiltonian as- 
sociated with a polarization P{t) induced by the external 
field £{t) is Hp{t) = —P{t) ■ £{t). Linear response the- 
ory now yields the Cartesian /-component of the induced 
inter-band polarization due to a weak external field £: 



Vi{t) 



dt'e{t-t'){[p{t'),pm)-m- (2) 



The retarded susceptibility tensor can be identified from 
Eq.(l): 



t'){[p„,{t'),p,m. 



(3) 



Following the standard line of attack in non-equilibrium 
theory [HJ], we first consider the causal (time-ordered) 
response function: 



-(T{P,„(t')P/(i)}) 



(4) 



where T is the time-ordering operator. In non- 
equilibrium, the causal response function is replaced by 
its contour ordered counterpart: Ximi^T^') ~* xFm('''j '''Oj 
where the complex-time variables r, r' reside on the 
Keldysh contour. Finally we obtain the retarded ten- 
sor by an analytic continuation using the Langreth rules 

We use Eq.(nl) to write the susceptibility as 



qk 



(5) 



X [{T,{al{r')b,{r')al{r)b,{r)}) 
+ {Tc{aliT')b,{r')bi{r)ak{T)}) 
+ {Tc{bliT')a,{r')al{r)b,{r)}) 
+ {Tc{bl{r')a,{r')bl{T)akir)})] . 

where Tc is the contour-ordering operator. In equilib- 
riuin, the two-particle correlation functions occurring in 
Eq.(||) would be found via the Bethe-Salpeter equation 
[p^ . In what follows, however, we shall consider the non- 
interacting limit of Eq.(^). This approach is motivated 
by the following considerations. The noninteracting limit 
will allow significant analytic progress, and the results, 
which we believe are interesting in their own right, will 
form the basis for any subsequent interacting theories. 
Secondly, experimentally it is known that in undoped 
semiconductor quantum wells excitons are quenched if 
the system is subject to intense FIR [^, thus implying 
that for the present situation Coulomb interactions may 
not play an equally dominant role as is the case in equi- 
librium situations. A quantitative assessment requires a 
non-equilibrium theory for the two-particle Green func- 
tions and will be a topic of our future work. For non- 
interacting particles we can use Wick's theorem to fac- 
torize the two-particle correlation functions. Thus, the 
non-equilibrium susceptibility can be expressed in terms 
of single-particle Green functions. The following Green 
functions are needed: 



gc{k,T;q,T') = -2(Tc{afe(T)at(r')}) 
9v{k,T;q,T') = 
gab{k,T;q,T') 



-^{TMr)bl{r')}) 
-^(T,{a,(r)6t(r')}) 



(6) 
(7) 
(8) 



and 
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9ba{k,T;q,T') 



-^(T,{6,(r)4(r')}). 



(9) 



We assume that the frequency, fi, of the FIR field is such 
that hn ^ eg. In typical experiments on III-V systems eg 
is of the order of eV, while is a few meV, so this con- 
dition is satisfied. Consequently, inter-band transitions 
due to the perturbing field can be ignored, and the Green 
functions related to Zener-effect, i.e., gab{k,T\q,T') and 
gbaik, t; q, r'), are neglected from this on. The first order 
non-equilibrium susceptibility reads thus 

xL(t,t') = -^^di{k)dm{q)[gcik,T;q,T')gy{q,T';k,T) 

q k 

(10) 



+9v{k,T;q,T')gc{q,T';k,T)]. 



The analytic continuation to real times is performed with 
the Langreth rules which state that if on contour 

C{T,T')^AiT,T')B{T',T), (11) 

then the retarded function on the real time axis is 

C"^(t,t') ^ A<{t,t')B''{t\t) + A'-{t,t')B<{t',t). (12) 
We thus have 

XlmM - -J5]di(fc)d™(fc)[5<(fc,t,0<(fc,t',t) 
k 

+g:{k, t, t')g<{k, t', t) + 5<(fc, t, t')gt{k, t' , t) 



with 



+glik,t,t')g<{k,t',t)], 

g<{t,t')=^{c^{t')c{t)), 
g^it,t')=zeit'-t){{cit),cHt')}}, 
g^{t,t') = -i9{t-t'){{c{t),cHt')}). 

We recall the following relations: 

[g<{tX}Y ^-9^{t',t) 
[g''{t,tT -5'-(i',i) 
[g'-{t,t')r ^g''{t',t) 



(13) 

(14) 
(15) 
(16) 



(17) 



For certain applications, e.g. Section III A, it is con- 



venient to introduce the center of mass variables T ~ 
{t' + t)/2 and t = t — t' |jl^. In terms of these variables 
the symmetry relations of the Green functions are: 



[g<(r,r)]* =-.g<(r,-r) 
[g%T,T)r ^g^iT,-r) 
[g'-{T,T)r ^g''{T,-T) 



(18) 



The retarded susceptibility expressed in center of mass 
coordinates is 



Y,Mk)dmik) 



x{[g<{k,T,T)g:{k,T,-T) 
+g<ik, T, T)gl{k, T, -r)] - h.c.}. (19) 



Note that i n equ ilibrium Ximi'^^''') — Ximi'^)- shown 
in Section [II A , the relevant quantity for continuous 
wave measurements at frequency uji is 



^i^XhniT,UJi) = Im 



dTe-'-xL(r,r) (20) 



to first order in fl/uji (here is the FIR frequency). Now, 
Xj„j(T, r) is a real quantity as is evident from (p^. Using 
the properties of the Fourier transform we obtain 

xUT,cui) = -^Y.diik)d„.{k) r ^{g<(fc,r,c.) 

x[g^{k,T,LU-coi)+gl{k,T,uj + coi)] 
+gv{k,T,uj) 

x[g-^ik,T,LO ^ LOi) + g:{k,T,uj + uJi)]}. 

(21) 

Since XimC^' ''') t^^sI, the imaginary part of its Fourier 
transform is obtained through 

imxL(r,^,) = hx\„^{T,uji) - xL(T, -c.,)]. (22) 

We can therefore write, in terms of the spectral functions 
A,{k, T, u:) = i[gl{k, T, oj) - g^^k, T, cu)] (23) 

and 

A,ik,T,u)^i[g:ik,T,to)-g'^ik,T,Lu)], (24) 

that 

ImxL(T,c.O = ^Y.di{k)dUk) ^{g<{k,T,^) 

LU + U!l)] 

+9vik,T,uj) 

x[A,{k,T,uj - uji) ~ A,{k,T,uj + uji)]y 

(25) 

The lesser functions can be expressed in the form 

5<(fc, T, to) = z/a(fc, T, L0)Aa{k, T, u) (26) 

where faik, T, lu) is a generalized particle distribution for 
particles of species a, and Aa{k,T,LL!) is the correspond- 
ing spectral function. In accordance with our assumption 
about no FIR field induced inter-band transitions, we can 
set /c(fc,T, w) = (zero occupation of conduction band), 
and that fv{k,T,uj) = 1 (all valence states are occupied). 
In the general case, e.g., when considering nonlinear ef- 
fects in the probing light field, one would have to find 
fa{k,T,uj) via, say, a Monte Carlo solution of semicon- 
ductor Bloch equations [|l^,^ or by a direct integration 
of quantum kinetic equations for g^{k,T,u;) With 
these assumptions the susceptibility reduces to 
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X {Ac{k, T,iu- Loi) - A,{k, T,iu + toi)} 
^^di{k)d,n{k) 



duj 
2^ 



Ayik,T,uj)A,ik,T,io + LOi). 

(27) 



The second equality comes about because we do not con- 
sider overlapping bands. Eq.(p7|), which is the central 
result of this section, expresses the fact that the non- 
equilibrium inter-band susceptibility function can be cal- 
culated from a joint spectral function, which is a convo- 
lution of the individual band spectral function. A sim- 
ilar result is known from high-field quantum transport 
theory there the field-dependent scattering rate is 
expressed as a joint spectral function for the initial and 
final states. 

In order to make a connection to the equilibrium case, 
we recall the exact identity g^^{k,uj) = iacq{k,uj)n(uj) 
{'n{u!) is the Fermi function), and obtain from Eq.(pq) 



ImxL(wi) 



-y 

k 



di{k)dm{k) 



'^\nc{uj)ac.cq_{k,Lo) 

X [a„,oq(fc, iO - UJl) - a„,eq(fc, iO + UJl)] 
x[ac,oq(fc,t^ - i^l) - ac,cqik,UJ + LUl)]y 

(28) 

Here nduj) is the conduction band electron occupation 
function, and nu(tt') is corresponding function for the va- 
lence band electrons. 



III. ABSORPTION COEFFICIENT IN TERMS OF 
THE TIME DEPENDENT DIELECTRIC 
SUSCEPTIBILITY 

The dielectric susceptibility x links the induced polar- 
ization V to the field £ via 



vit) 



dt'xit, t')E{t' 



The wave equation for light is then 



1 d'^V{t) 







(29) 



(30) 



where P(i) — £{t) + 47r7'(t). The absorption coefficient 
a{uj) is defined as the inverse of the length which light 



has to traverse in the medium at frequency uj in order for 
the intensity of the light to decrease by a factor of 1/e. 
In equilibrium X'(a') = [1 -I- 47rx(a;)]f (w) = t{uj)£{uj) and 
the absorption coefficient p8| becomes 



a{ijj) — Attuj 



Inix(^) 
cn{uj) 



(31) 



Here n^(w) = i[Ree(a;) -|- |e(a;)|] is the refraction coeffi- 
cient which usually depends only weakly on u. In non- 
equilibrium this analysis must be generalized, and we 
consider two special cases: (i) monochromatic continu- 
ous wave measurements, and (ii) white light short pulse 
measurements. 



A. Continuous wave measurements 

Consider a system out of equilibrium which is probed 
by a light field (which is assumed to be weak) of frequency 



£{r, t) = £o exp[i(r • k — uit)]. 



(32) 



The polarization can then be expressed as 

V{t) = £{r,t) diV"'(*"*')x''(t,i') • (33) 



This form is suggestive: it is advantageous to express 
in terms of the center-of-mass and difference coordi- 
nates, X^it^t') ~^ X^i^i^ + t'),t — t'). The characteristic 
time-scale for the center-of-mass time is set by the "slow" 
frequency f2, while the difference-time varies on the scale 
of the "fast" frequency cui. We thus gradient-expand: 

x-{^it+t'),t-t')^^ rit,t-t')+^it' -t)r' {t,t-t')+- ■ ■, 

where the prime indicates differentiation with respect to 
the slow temporal variable. Substitution in ( p3| ) then 
yields (we introduce a new variable t = t — t') 

V{t)^£{r,t) J dTe^'^'^[r{t,T) + i-^T)r'{t,T) + ---] 
= £(r,i)[x'-(t,^0 + ^|^r(i,^0 + ---] 



£{r, t) exp 



i d 



,2 n 



2 dtdioi 



(34) 



Eq.(p4[) can now be used in the Maxwell equation; note 
however that upon Fourier-transforming the dominant 
frequency comes from £{t) and we can keep t in x(i,a;;) 
fixed. The slow time-variation will from this on be indi- 
cated by T. Proceeding as in deriving the static result 
(^), we identify the time- dependent absorption coeffi- 
cient 



cnxyuj) 



(35) 
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If the driving force is periodic in T (the harmonic time- 
dependence due to a FEL-laser is an important special 
case), then the average absorption coefhcient is 



a{uj) 



1 



period J period 
1 



AT aT{Lo) 



(36) 



dTAnuj 



period J period 



Imx'"(T, Uj) 
cnT{uj) 



to all orders in il/uj. We stress that here x^{T,uj) is 
Fourier transformed with respect to the difference vari- 
able T. Below we shall represent numerical examples for 
the generalized absorption coefficient. 



B. Short white light pulse measurements 

Consider now an instantaneous measurement per- 
formed on a non-equilibrium system: at some specific 
time t = tm the system is probed with a weak pulse 
whose duration is short compared to the characteristic 
dynamics of the system. We approximate the pulse with 
delta function in time: 



(37) 



In principle, construction of such a pulse would take in- 
finite energy due to its time dependence. The pulse is 
therefore hardly "weak" . When we refer to the pulse as 
weak we assume that Eq ^ 1, that is the intensity of 
the light is small at all frequencies. Using (|3^) in the 
Maxwell equation yields dispersion relation 



ui 



[l + 47rx'^(a;,i™)]. 



(38) 



This dispersion relation looks quite similar to the one ob- 
tained in the previous section. The difference is that here 
X^iuJ,tm) is Fourier transformed with respect to t' , not 
the difference variable r. In the present case we obtain 
the time dependent absorption coefficient 



at{w) = Attuj 



cnt{u)) 



(39) 



For examples of experiments which probe systems in this 
manner see, e.g., Refs. |9|JlO|]. 



C. Differential Transmission Spectrum 

Consider a sample of thickness L; then the ratio of 
the intensity transmitted through the sample with its 
initial intensity is T{u)) — exp[— Q:(ti;)L] where a{uj) is 
the absorption coefficient of the sample. Experimental 
setups for measuring the change in absorption due to 
externally controlled perturbations commonly measure 
the differential transmission spectrum (DTS) defined by 



DTS{uj) 



Tjoj) - Tpjuj) 
To{lo) 



(40) 



Here T(w) is the transmission with the perturbation 
present and Tq{uj) is the transmission through the un- 
perturbed sample. Below we give examples of DTS{uj) 
in non-equilibrium situations. 



IV. SINGLE PARTICLE GREEN FUNCTIONS 
AND THE SPECTRAL FUNCTIONS 

In this section we determine the single-particle Green 
functions and their associated spectral functions. We 
show that, under conditions specified below, that the con- 
volutions of the spectral functions, encountered in Sec- 
tion |l|, result in effective single band spectral functions. 



A. The single particle Green functions 

Let A be the vector potential which defines the FIR 
field. Considering harmonic, translationally invariant ex- 
ternal fields we choose 



A{t) = -E 



(41) 



which represents the physical uniform electric field 
Ecos{ftt). The Hamiltonian describing the two bands 
is then 

H = Y.{ec[k- eA{t)]alaj: + e^k - e^t)]^^^} . (42) 
fe 

The Dyson equation for the retarded/advanced free- 
particle Green function is 

(^ih-^^~e^[k-eAit)]^g:/^ik,t,t')=Sit-t') , (43) 

where a G {c, v} is the band index. This equation is 
readily integrated with the solutions 



fJ^ik,t,t')^T^9{±tTt')e^p{ 



dsea[k — eA(s)]|, 
(44) 

and the spectral function A ~ i{g''' — g"') becomes 



A. 



,{k,t,t') = iexp{ - dse^ik - eA{s)]y (45) 
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B. Convolution of the spectral functions 

According to Section |^ the susceptibility is obtained 
through the trace of a convolution of the spectral func- 
tions. We shall now show that within the present model 
the convolution of spectral functions results in a new ef- 
fective single-band spectral function. 

In terms of the center of mass variables, t = t — t' and 
T = {t + t')/2, we write the spectral functions as 



A^{k,T,uj) = i / dre*' 



(46) 



X exp 



T-t/2 



Then the convolution 



—Ae{k,T,uj)A^{k,T,u - Loi), (47) 



of the two spectral functions becomes 

h{k,T,uji) = - dre^'^'^exp] -- / As (48) 

x{e,[k-eA{s)]-e4k-eA{s)])y 
In the case of parabolic bands we have 

h^k^ 



n^k'^ 

Zirie 



2m/i 



(49) 



Here me is the electron mass, to^ is the positive hole 
mass and Cg is the band gap. We define a single effective 
band for the system: 

- - - ?j2^2 

eeff[fc] EE ec[k] - e„[fc] = + eg, (50) 

2mcff 

where rrieff = memh/{nie + nih) is the effective reduced 
mass. Thus, the effective band is parabolic like the orig- 
inal bands but with their reduced mass. It is there- 
fore evident that the convolution, writing Acs{k, T, lui) = 
b{k,T,uji), is a spectral function for a parabolic band. 



A,s{k,T,uj) = ^ I dre" 



exp{-- 



(51) 



T+t/2 



T~t/2 



dsecs[k ~ eA(s)]}. 



In the case of tight-binding minibands for a type I super- 
lattice (with period a) we write the bands as 

1 f?k^ 
ec[fc] = -^\c cos(afc|| ) -I- ^— i (52) 

1 h^k^ 
(■v{k\ = -^\v'^Q^{ak\\) - f-g, (53) 



where Ac is the electron miniband width, is the cor- 
responding bandwidth for the holes; fc|| is the (crystal) 
momentum component parallel to the growth direction of 
the superlattice and k\^ is the magnitude of the compo- 
nent perpendicular to the growth direction. The effective 
band thus becomes 



EcffW = ^Aeff cos(aA;||) 



2 1,2 



rfk 



2top 



(54) 



where Acff = Ac + A^, is the effective bandwidth, which 
again is of the same form as the original bands. This 
shows that also for superlattices the convolution ( |47| ) 
leads to an effective spectral function of the original form. 

In terms of the effective spectral function the imagi- 
nary part of the susceptibility can be written as 



did„ 
~2h 



^^cff(fc, T, uji), 



(55) 



where we assume that the dipole matrix elements are 
/c- independent. In equilibrium the trace of the spectral 
function yields the density of states for the system. Anal- 
ogously, the generalized time- dependent density of states 
(GDOS) [| is defined as 

p{T,Lu,)^-Y,AMk,T,L,i), (56) 



allowing us to write the absorption coefficient as 

2TT^UJl\d\^ 

aT[u}i)Ki — p[T,u;i). (57) 

cnn 

For the remainder of this work we shall investigate the 
properties of p(T,uji) for various systems. 



C. Gauge invar iance 

To conclude this section we briefly comment on gauge 
invariance. We might have, from the outset, chosen 
to work within the gauge invariant formulation which 
has been developed in of high- field transport ||2^,^,|l2| . 
Considering translationally invariant systems, correlation 
functions are made gauge invariant with the transforma- 
tion 



k + 



t-t' 



dsA{s). 



(58) 



However, the absorption coefficient follows from a trace 
operation (^6|), which makes the transformation ( |5^ ) re- 
dundant: a simple change of variables when performing 
the trace undoes (^8|) and proves that our formulation of 
the absorption is gauge invariant. 
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V. PARABOLIC BANDS 

In this section we shall investigate the properties of the 
generalized density of states for systems which can be 
effectively described by Hamiltonians yielding parabolic 
bands, be it in one, two or three dimensions. We write 
the effective single-band dispersion as 



= H Co 



2moff 



(59) 



For convenience we write m = nieS and set eg = which 
shifts the energy axis such that the reference point is the 
band gap energy. We calculate the generalized density of 
states from 



p"^(r,e) 



dre"^/'"V"-°(T,T), 



(60) 



where 



I - - / dse[k-eA{s)]}. 

(61) 



With the vector potential ( [4l| ) one obtains explicitly that 
p"^(T,r) = -y ^exp{-z[(e, + e,)r/n 



ek ■ E ^ 



n 



cos(217r) sin(r2T)]|. 



(62) 



Here — h^k^ /2m and we have defined the fundamental 
energy scale 



— hujt — 



(63) 



The energy e f can be interpreted classically in the follow- 
ing way: consider a classical particle with charge e and 
mass TO subjected to an electric field E{t) = Ecos{nt). 
From Newton's equation of motion one finds that the 
mean kinetic energy of such a particle equals e/. 

In order to perform the Fourier-transform (pO) we uti- 
lize the identity pi] 



exp(ia; sin 9) = Jn{x) exp(m0). 



(64) 



where J„(a;) are Bessel functions; we shall henceforth 
write J2n = '^nL-oo ^'^ simplify the notation. The gen- 
eralized density of states becomes 



xJi+,(^cos(2rir) 



The dimensionality is entirely contained in the remain- 
ing momentum integration J d"fc/(27r)"~-'^. We note that 
Eq.(|65|) imphes a shift of the absorption edge by e/. The 
term Ihn in the Dirac-delta function gives rise to pho- 
tonic side bands. Since J2i{x) is an even function, the 
density of states is invariant under the transformation 
E —* —E, as expected. In the following subsections we 
shall consider the ID, 2D and 3D systems separately and 
show how the density of states smoothly evolves from a 
low field intensity regime into a high field intensity regime 
making the non-linear effects of the THz field apparent. 



A. Generalized density of states, ID 

The density of states for a single-mode ID-system 
( "quantum wire" ) in the absence of external fields is 



1 /2to 



1/2 



-1/2 



9{e). 



(66) 



In the presence of a external strong oscillating field we 
get from (^ that the GDOS is 

/OO 
dkd[e-ek-ef + inn] 

( \/32e/efc . 



''''' \-Ml 



■ sin(17r) 



xj.,,(i^cos(2i7r)^ 

= ^ r]^ (T, e-ef + lhn)pl^ (e - £/ + Ihn) , 

(67) 

where the side-band weights are 

r-(T,e)=^J,,Y^^sin(17T)j 



XJ/+J (^cos(2r!r) 
We note that in the limit e/ ^ 



(68) 



(69) 



and p^^{T, e) /oj^(e), as expected. If fl is in the THz 
regime then most experiments would probe the time av- 
eraged absorption. The time-average of the side-band 
weights is calculated from 



:;1D 



2^ ds 



'■2-3 



sin(s) 



(70) 



(65) 

which yields for I odd: 
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E 



X J2j+i+i (^cos(2s) 



(71) 



and for I even 



At the onset of side band I the side-band weight is 



(0) 







if / odd 



J;^2 / /^^) if ' even, 
where we used the identity p3] 



2Tr 



Aej2i{acos0) = 2Tijf{a). 



(72) 



(73) 



(74) 



This shows that processes involving an odd number of 
photons of the THz field are strongly suppressed. In Fig. 
|l| we illustrate Pave(^) ^ range of values of ef/hfl. 
In the figures we write Ce = hil. We observe all the 
signatures of the Dynamical Franz-Keldysh effect (DFK) 
1^: Stark- like blue shift of the main absorption-edge by 
e/, formation of side-bands at + ± Nhfl, and finite 
absorption within the band-gap. 



B. Generalized density of states, 2D 

Several authors have considered fields perpendicular to 
the quantum well, cf. | |25[ | and references in these papers; 
here we focus on the situation where the electric field is 
in the plane of the two-dimensional electron gas. In such 
a system with no external field the density of states is 
constant. 



(75) 



With a harmonically oscillating field we obtain from (p3) 



where the side-band weights are 



„2D 



(78) 



xJi+,(|^cos(217T)). 
Identical arguments as in the ID case lead to 



(0) 







if I odd 



Jfl^i^f/hil) if / even. 



(79) 



i.e., the same result as in the ID-case. 

As in the ID case we have numerically investigated the 
time averaged GDOS plg{e) = £ j^""'^ dT p^^{T,e). In 
Fig. 1^ we illustrate Pa™(£) fo'" various values of ef/hfl. 
Again, as in the ID case, we observe all the characteristics 
of the DFK B. Finally, Fig. | shows the differential 
transmission (DTS) signal. 



C. Generalized density of states, 3D 

Absorption in bulk semiconductors subject to THz 
radiation was considered already long time ago by Ya- 
coby He studied transition rates between bands by 
investigating approximate solutions to the correspond- 
ing time-dependent Schrodinger equation. He concluded 
that transitions occur in the gap and noted reduced rates 
above the gap, both in agreement with the present work. 
General quantitative results, were not presented. The 3D 
field-free density of states is 



„3D, 



With the external field the density of states becomes 
p3D(T, e) = ^ rf^{T, e-ef+ inVL)pl°{e - e/ + IhVL) , 



(81) 



with the side-band weights 



,-(T,.).|:y^ dkkj^ -S[e-e,-e, + inn] .-(^^ ^) ^ ^ eos(2^)T)) 



xJ2j ( ^^^^cos6lsin(f]T) 



xJi+j (|^cos(2r!r) 



x/ d?7J2j I -^4^sin(J]T)r;] . (82) 







(76) 



The integrals in (|7^) are again performed using ( |7^ ) and 
writing the result in the side-band picture we obtain ||2^ ] 

p2^(T, e)^Y. ^'''(^' '-'f+ ^mpVi^ -^1+ m 
I 

(77) 



Again, we have 



if I odd 

Jfl^{€f/h^) if / even. 



(83) 



In Fig. ^ we illustrate Pave(£) fo'' various values oief /fiVL\ 
the DTS signal for the 3D case is shown in Fig. ||. 
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D. Summary 

The main physical consequences of the THz field on 
linear absorption spectrum for systems with parabolic 
dispersion can be summarized as follows. The dynamical 
modifications of the absorption spectrum (i) Appear near 
the absorption edge; (ii) They extend a few Ce = hVL 
around the edge, and (iii) They are most pronounced 
when LOf /D, is of order unity. 

If H, is in the THz regime, and fields like those attain- 
able with free electron lasers are considered [||j2^], then 
ujf/fl « 1 and the fine structure extends over an area of 
several meV. Consequently, an experimental verification 
of these effects should be possible. 



VI. SUPERLATTICES 

According to the semiclassical Bloch-Boltzmann the- 
ory of transport, a uniform electric field causes charge 
carriers in a periodic potential to execute a time-periodic 
motion with frequency ujb = eaE/h, where a is the lat- 
tice periodicity. Conditions for observing these Bloch os- 
cillations are much more favorable in superlattices than 
in ordinary bulk materials, and recent years have wit- 
nessed an intense research effort culminating in the ob- 
servation of Bloch oscillations j2^. In ac-fields a phe- 
nomenon called dynamical localization may occur: if the 
parameter 7 = aeE u IMl equals a zero of Jq, the average 
velocity vanishes [g9| . In this section we investigate how 
dynamical localization [^||29|-|3^ manifests itself in the 
free particle absorption spectra. Recently, Meier et al. 
p3| presented results of a detailed numerical solution of 
the semiconductor Bloch equations, including excitonic 
effects, and found that at dynamical localization the rel- 
ative motion exciton wave function changes from a 3D- 
character (i.e., localized in fc^-space) to a 2D-structure 
(extended in fcz-space) , and below we shall illustrate how 
the same phenomenon reflects itself in the present ana- 
lytic study of free-particle properties. 



A. Generalized Density of States 

The starting point for our analy sis is the effective dis- 
persion ( p^ introduced in section IV B which we repro- 
duce here for the convenience of the reader, 



ecff[fc] = ^Acff cos(afc||) 



2mes 



^9- 



(84) 



Henceforth we put e,, = and drop the "eff" sub- 
script. We consider the effect of the THz field described 
by the vector potential A{t) = —Esm{nt)/n where 
E = (0,0, In accordance with Section IV we cal- 

culate the generalized density of states from 



p^'(T,r) 



xexp[/(r,T)] , 



00 j a 

dei / dfciie-"^^/'' 



(85) 



I{T,T)^~i— dscos[aA:|| -|-7sin(r2s)]. (86) 

2h Jt-t/2 

We evaluate the integral within the exponential using the 
identity (Q) with the result 

^^^'^^ " ^{ cos(afc|| )[C(f7r) + ^0(7)^] 

+ sin(a/c||)5(17r)| , (87) 



where 



rf \ o "^2n(7)cos[2nl7r] . 

C[z)~2y sm[n2;J (88) 

Tl 

n—1 

,_^^ J2n-ih)sm[{2n-l)nT] 



n=l 

X sin[(n — 1/2)2;] 



(89) 



We have suppressed the explicit dependence of fJT and 7 
in C{z) and S{z) for simplicity. Note that both of these 
functions are anti-symmetric in z, i.e., C{~z) — ~C{z) 
and S{-z) = -S{z). The identity [||| 

/ — exp{mcos6'-^ j6sin6'} = Jo{Va'^ + b^) (90) 
Jo 27r 

is the key to the next step in the evaluation of ( p5[ ) and 
allows us to write 



Jo 

where we have defined the kernel 



de_L e-"^^/'^/C(17T) (91) 



m = Jo V(C(z) + Ml)zr + S^z)^ . (92) 

Also here we have suppressed the explicit dependence of 
flT and 7. In distribution sense we can write 



de^e-"^^/'* = 



T-iO+' 



(93) 



where 0+ indicates a positive infinitesimal. This expres- 
sion allows us compute the Fourier transform of (|9l|): 



2TT^h'^C 



, sin(eT/?i) _ ,^ , 
dr — ^ ' ' 1C{Q.t) + TT 



(94) 



In what follows we shall examine several properties of 
this result. 
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B. The field- free limit 



D. Dynamic localization 



In the limit of vanishing field strength we have 
lim/C(.) = Jo(^ 



(95) 



Using the identity 



X 



[(3 > 1] 

sm{(3x)Jo{x) ^ { 2arcsin/3 [l3^ < 1] (96) 

[f3 < -1] 



we obtain the density of states for a tight-binding super- 
lattice, 



irh a 



1 [e > A/2] 

i arcsin(2e/A) + 1/2 [|e| < A/2] (97) 

d [e < -A/2] 



which is familiar. 

C. The static limit 

In the limit ^ one obtains 

lim JCiyiT) = Jo ( — ^ sin(tJBT) 
n^o \2nujB 

Recalling the identities ||2^ 

Jo (z sin a) = ^ jJ(z/2) cos(2fcQ;) 
j 

and 



(98) 



(99) 



dx cos(ax) sm{Px) /x — 7T[6{a + P) — 0{a — f3)] , 

(100) 

we obtain the density of states 



(101) 



This result coincides with the one obtained in Refs. 
p^ , ^ , which study both theoretically and experimen- 
tally the effects of strong static fields on the absorption 
in super lattice structures. They conclude that the step- 
like behavior of ( |101| ) is due localization in the growth 
direction (Wannier-Stark localization). 



As seen in the previous subsection the signature of lo- 
calization in the growth direction in a superlattice is a 
step-like behavior of the density of states. This is in- 
tuitively clear since the density of states for a 2D sys- 
tem ( [TSl ) is constant. We therefore expect the density of 
states to be composed of a step function for each well 
the states extend into, with weight relative to the occu- 
pation in that particular well. We shall now show that if 
"^0(7) = 0, i.e., the conditions for dynamical localization 
are met, then the GDOS indeed is of this kind. The ar- 
gument runs as follows. If Jo (7) = then the kernel ( |9^ ) 
is periodic in z with period 27r. Furthermore, the kernel 
is an even function: JC{z) = lC[—z). We can therefore 
formally write 



/C(f)T) =^/CjCOs(fcnr), 



which is of the same functional form as in the static limit, 
Eg. (p8|) . Consequently, we may conclude that the gener- 
alized density of states must be of the form 



dyn.loc 



(e) 



jhn), 



(102) 



i.e., it is a superposition of step- functions. The weights 
Kj, however, must be evaluated numerically, and exam- 
ples are given in the next section. It is important to note 
that the "step-length" in the ac-case is determined by 
the frequency of the THz-field in contrast to the static 
case, where it is determined by the field strength. The 
field strength enters the density of states only through 
the weight-fac tors Kj. 

The result (102) suggests that it should be possible to 



probe dynamic localization by photo-absorption: when 
the appropriate conditions are approached, the absorp- 
tion coefficient should change qualitatively from a generic 
smooth behavior to a sharply defined step-like structure. 
The number of distinct steps appearing in the spectrum 
is determined by the ratio X/Ml which is also a measure 
of the number of wells the localized states span. This 
is fully consistent with the results of who consid- 
ered a miniband of width A = 21meV and FIR frequency 
TiU, = 20meV, which corresponds to a single step, and 
hence maximum binding energy of the corresponding ex- 
citon which would be mostly confined to a single well. 



E. Numerical results 



We again focus our numerical study to the time 
averaged generalized density of states Pavo(^) — 
n_ j2^/n ^j,^si(j^ pigg^ I j^^j I ^Yiow the ab- 

sorption spectra for a superlattice with effective band- 
width A = 3AhQ. The numerical results confirm the 
expectations of the previous section: when the 7 = 
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aei5||/?if2 approaches the first zero of Jq, which occurs 
at the argument value of 2.4048.., the gradually evolving 
replicas of the zero-field density of states converge into 
plateaus of finite width. The exactness of the plateaus 
can be judged from Fig. 0: at 7 = 2.4048.. the Hne 
joining the the steps appears near vertical. Finally, in 
Fig. H we show the DTS spectra at dynamical localization 
(DL) and non-DL conditions. There are two character- 
istic differences: (i) Outside the zero-field miniband DL 
leads to a step-like structure in contrast to the smooth 
behavior found otherwise, and (ii) Inside the miniband 
the DL-spcctrum distinguishes itself by its sharp jagged 
structure. 



VII. CONCLUSIONS 

We have presented a theoretical formulation of lin- 
ear photo-absorption for samples under strongly non- 
equilibrium conditions. Typical non-equilibrium agents 
would be THz-radiation from free-electron lasers, or 
ultra-short pulse measurements of transient effects. In 
the present work noninteracting carriers are considered, 
but the formulation allows an extension to Coulomb in- 
teractions, which will be addressed in our future work. 
Two central concepts emerge from our analysis: a gen- 
eralization of the density of states into time-dependent 
conditions [GDOS defined in Eq.(|56|)], and photonic side- 
bands, which form a convenient framework for discussing 
the various features of the absorption spectra. 
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FIG. 1. Time averaged generalized density of states 
for ID-system shown for a range of FIR-intensities, 
€f/Ml G [0,2]. The band edge and the side-bands display 
a blue-shift, which scales linearly with the intensity. Absorp- 
tion extends below the bandgap (dynamical Franz-Keldysh 
effect). 




(e-eg)/e, 

FIG. 2. The time averaged GDOS for a 2D-system for 
a range 
of FIR-intensities, e; /TiVi = (0.2, 0.5, 0.8, 1.1, 1.4, 1.7, 2.0). At 
low intensities one observes a Stark-like blue-shift of the band 
edge as well as finite absorption within the band gap. With 
increasing intensity side bands emerge at e = eg + e/ ± 2TiQ.. 
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-20246 
(e-eg)/ee 

FIG. 3. The DTS signal for a 2D-systcm for a range of 
intensities, tf/hQ. = (0.2, 0.5, 0.8, 1.1, 1.4, 1.7, 2.0). 




-2-1012345 



FIG. 4. The time averaged generalized density of 
states for a 3D-system for a range of FIR-intensities, 
tf/hQ, = (0.0, 0.2, 0.5, 0.8, 1.1, 1.4, 1.7, 2.0). 




FIG. 5. DTS-spectrum for a 3D-system. 



-3 -2 -1 1 2 3 
(e-eg)/ee 

FIG. 6. The time-averaged GDOS for a superlattice with 
effective band-width A = 3.4?!,f2. Dots: 7 = 1.5; dashes: 
7 = 2.4. The proximity of dynamical localization occurring 
at 7 = 2.4048... reflects itself in the step- wise structure of the 
dashed curve. 




FIG. 7. The time averaged GDOS for a superlattice with 

A = 3.4?!,f2 as a function of FIR-intensity. At low 7 side bands 
are observed, which merge at 7 = 2.4048.. corresponding to 
dynamical localization. 
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(e-eg)/e, 

FIG. 8. The differential transmission signal for super- 
lattice structure with effective mini bandwidth A = 8h^l. 
s Solid line: 7 = 2.4; dashes: 7 = 2.0. Outside the 
zero-field mini-band a step-like behavior is seen, while in- 
side the mini-band DTS for dynamical localization develops 
a jagged shape in contrast to the smooth behavior for the 
extended state. 
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